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A search for the exclusive radiative decays B — > p(770)7 and B° — > o;(782)7 is performed on a 
sample of about 84 million BB events collected by the BABAR detector at the PEP-II asymmetric- 
energy e + e~ storage ring. No significant signal is seen in any of the channels. We set upper 
limits on the branching fractions B of B(B° -» p°j) < 1.2 x 10~ 6 , B(B + -> p+j) < 2.1 x 10~ 6 , 
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and B(B° -> ury) < 1.0 x at 90% confidence level (C.L.). Using the assumption that T(B -> 

fry) = F(B + -» p+7) = 2 x r(B° -» p°7), we find the combined limit B(B -> P7) < 1.9 x 10" 6 , 
corresponding to B{B pj)/B{B K*j) < 0.047 at 90% C.L. 

PACS numbers: 12.15.Hh, 11.30.Er, 13.25.Hw 



Within the Standard Model (SM), the decays B -> /ry 
and B° — > ujj proceed primarily through an underlying 
b — > c?7 electromagnetic "penguin" diagram that con- 
tains a top quark in the loop These processes are 
analogous to the B — > X*7 process mediated by the 
b — > 57 transition, but with the final-state s-quark re- 
placed by a d-quark, and the relevant element of the 
CKM matrix changed from V ts to Vtd- There may also 
be contributions resulting from physics beyond the SM, 
such as supersymmetry [2|. Recent calculations of the 
branching fraction in the SM indicate a range B(B + — » 
p+7) = (0.9 - 1.5) x 10" 6 0,0. The range is due both 
to uncertainties in the value of Vtd and to uncertainties 
in the calculation of the relevant hadronic form factors. 
The rates for B° -> p°7, B+ -> p+7, and B° -> loj 
are related by the quark model, such that we expect 
T(B+ -» p+ 7 ) w 2 x r(B° -» p° 7 ) w 2 x r(B° -» w 7 ). 
Previous searches have found no evidence for these 
decays, nor any other b — > c?7 processes. 

The analysis uses data collected by the BABAR detector 
6| at the PEP-II asymmetric-energy e+e~ storage ring 
j|. The data sample consists of 84.4 ± 0.9 million BB 
events corresponding to 78 fb -1 on the T(4S) resonance 
( "on-resonance" ) , and 9.6 fb -1 recorded 40MeV below 
the T(AS) resonance ("off-resonance"). 

The BABAR detector consists of five subdetectors. 
Charged particle trajectories are measured by a combi- 
nation of a five-layer silicon vertex tracker (SVT) and 
a 40-layer drift chamber (DCH) in a 1.5-T solenoidal 
magnetic field. Photons and electrons are detected in a 
CsI(Tl) electromagnetic calorimeter (EMC), with photon 
energy resolution o E jE = 0.023(£/ GeV)" 1 / 4 © 0.019. 
A ring-imaging Cherenkov detector (DIRC) is used for 
charged-particle identification. The magnetic flux return 
is instrumented with resistive plate chambers to identify 
muons. 

The decay B — > pj is reconstructed with p° — > 7r+7r~ 
and p+ — > 7T+7T , while B a — > ui^f is reconstructed with 
lo — > 7r+7r _ 7r°. Charge-conjugate channels are implied 
throughout this paper. Background high-energy pho- 
tons are produced primarily in continuum u, d, s, and 
c quark- antiquark events through ir° fr\ — ► 77 decays or 
via initial-state radiation. The reconstruction uses quan- 
tities both in the laboratory and T(45) center-of-mass 
frames, where the latter are denoted by an asterisk. 

The primary photon in the B decay is identified as an 
energy deposition in the EMC. The deposition must meet 
a number of criteria (described in detail in our paper ||| 
on B — > K*j) that are designed to eliminate background 
from charged particles, hadronic showers, and ir° and 77 



decays. 

As in Ref. , the charged tracks used in identifying the 
p/uj meson are well-measured tracks with a momentum 
transverse to the beam direction greater than 0.1 GeV/c. 
A charged pion selection based on dE/dx measurements 
in the SVT and DCH, and on Cherenkov photons recon- 
structed in the DIRC is used to reduce backgrounds from 
the b —> S7 processes by rejecting charged kaons (e.g. K + 
from B° —> K*°j). Figure^a) shows the particle identi- 
fication performance measured with a control sample of 
D*+^ £>°(-> K-tt+)tt + decays. 

Neutral pion candidates are identified using two pho- 
ton candidates reconstructed in the calorimeter, each 
with energy greater than 50 MeV. The invariant mass of 
the pair is required to satisfy 115 < m 77 < 150 MeV/ c 2 , 
which removes pairs whose invariant mass differs from 
the true m^o by more than about 3 times the experimen- 
tal resolution. A kinematic fit with m 77 constrained to 
m^o is used to improve the momentum resolution. 

A p° candidate is reconstructed by selecting two iden- 
tified pions that have opposite charge and a common 
vertex. The p+ candidates are obtained by pairing ir° 
candidates with an identified charged pion. We select p 
candidates with invariant mass within 250 MeV/c 2 
of 771 p = 770MeV/c 2 @ and momentum 2.3 < p* w < 
2.85GeV/c. The lo candidates are reconstructed from 
combinations of oppositely charged identified pions with 
a common vertex and 7r° candidates with invariant mass 
m I+rT o within 23MeV/c 2 of to w = 783MeV/c 2 and 
momentum 2.4 < p* +7r _ 7r0 < 2.8GeV/c. The m^+^-^o 
resolution is slightly poorer in data than in Monte Carlo 
(MC) simulation. The resulting change in signal effi- 
ciency of the rriir+Ti-Tr selection is accounted for as a 
systematic error in the signal efficiency. 

The photon and p/ui meson candidates are combined 
to form the B meson candidates. We define AE* = 
Eg - E£ cam , where £* eam is the energy of each beam 
and E* B = E* + E*, is the energy of the B meson can- 
didate. The signal candidates are centered at AE* = 
with resolution of about 50 MeV and a tail towards neg- 
ative AE* due to the asymmetric energy response of the 
EMC. We also define the beam-energy-substituted mass 

toes = \/ ^bcam — p'b ' wnere p'jg is the momentum of 
the B candidate modified by scaling the photon energy to 
make E*+E*, u - E£ eam = 0. This procedure reduces the 
tail in the signal toes distribution that results from the 
asymmetric calorimeter response. The signal candidates 
peak at toes = m B with a resolution of about 3 MeV/ c 2 , 
dominated by the beam-energy spread. 



FIG. 1: a) Pion efficiency and kaon misidentification rate of 
the charged pion seiection for 0.4 GeV/c < pi ab < 4.0 GeV/c, 
where pi a t is the track momentum in the laboratory frame, 
b) B° — > p°7 neural-network output for MC-simulated events 
with comparison to data control samples. Events with neural- 
network output greater than 0.9 are selected, as indicated. 



We consider candidates in the "fit region" —0.3 < 
AE* < 0.3 GeV and 5.20 < to ES < 5.29GeV/c 2 . For the 
small fraction of events (8% for MC B° — > p°j events) 
in which more than one B meson candidate satisfies all 
the cuts, the candidate with the smallest value of |Ai?*| 
is selected. 

We construct a number of variables that distinguish 
the signal from the continuum qq background. As 
in Ref. 8], we calculate the thrust angle 8j>, the B- 
production angle B , and the helicity angle Oh- For 
B° — > W7, 9h is defined as the angle between the normal 
to the decay plane of the ui and the flight direction of the 
B meson, both computed in the to rest frame. We also 
calculate several additional discriminating variables. The 
energy flow of the event excluding the B-meson daugh- 
ters in 10° cones centered on the photon-candidate mo- 
mentum provides discrimination between the jet-like con- 
tinuum background and the more spherical signal events. 
For suppression of the initial-state radiation background, 
we consider R' 2 , the ratio of second- to zeroth-order Fox- 
Wolfram moments jl0| in the frame recoiling from the 
photon momentum. We define the net flavor content as 
Y^i \N* — AT - 1, where are the number of e^, /Lt ± , K^, 
and slow pions of each sign identified in the event [llj |. 
On average, BB events have larger net flavor than con- 
tinuum events. In the B — > p 0r y and B° — > ujj analyses, 
we use the seperation along the beam axis of the B-meson 
candidate vertex and that of the rest of the event. Due 
to the finite B lifetime, this should be larger in magni- 
tude in BB events than in continuum background. In 
the B° — ► oj'y analysis, we use the uj Dalitz angle Qd, 
which is defined as the angle between the 7r° and the tt + 
in the 7r + 7r~ rest frame [l.2j; cos6*£> follows a sin 2 0d dis- 
tribution for true u> decays, as opposed to the uniform 
distribution of combinatorial background. 

The background-suppression variables are combined 
into one discriminating variable via a neural network, 
which responds non-linearly to the input variables and 
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FIG. 2: AE* vs. m E s fit regions for a) B° p°-y, b) B + -» 
p + 7, and c) B° — > uj'y candidates. The boxes indicate the 
regions where signal events would appear: —0.2 < AE* < 
0.1 GeV and 5.27 < m E s < 5.29 GeV/ c 2 . Assuming B{B° -> 
p° 7 ) = -> p+7) = B(B° -> w 7 ) = 10" B , we expect 

9.9, 12.1 and 3.4 signal events in these regions, respectively. 



exploits correlations between the variables |13|]. A sepa- 
rate neural network is trained for each mode. 

The output for the neural network trained for £>° — > 
p°7 is shown in Fig. Q^b), where the MC simulation 
of the continuum background is compared with the off- 
resonance data, and the output for MC-simulated B° — > 
D~ir + decays is compared with B° —> D~ir + decays re- 
constructed in the on-resonance data. The latter com- 
parison provides a cross-check of those input variables 
that depend on the properties of the other B meson in 
the event. This includes all of the variables except for 
Oh and Op, which, for this check, are modeled using the 
signal MC distributions. 

To suppress the continuum background, we make a 
selection on the neural-network output that is optimized 
for minimum statistical error as determined using MC 
samples of signal and background. The efficiency of this 
selection for the B — > Dir control sample differs slightly 
between the data and MC. We account for this difference 
as a systematic error in the signal efficiency. For B + — > 
p + 7, we also require | cos Oh \ < 0.6 to reject B + — > p + 7T° 
events, which have a cos 2 Oh distribution, as opposed to 
the expected sin 2 Oh distribution of the signal process. 

After applying the neural- network, cos Oh, and fit- 
region selection to the on-resonance data, 449 events re- 
main in the B° — > p°7 data, 480 events for B + — > p + 7 
and 54 events for B° — > wy. MC studies indicate that 
about 90% of the background in these samples comes 
from continuum events, and only about 10% from BB. 
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Mode 


Yield 
(Events) 


Bias 
(Events) 


Upper Lim. 
(Events) 


e B 

(%) (io- 6 ) 




P°7 4.81^ 


[-0.5,0.8] 


12.4 


12.3 0A±%% 


B+ -> 




[-0.1,2.0] 


15.4 


9.2 0.7l° ;l 


B° — » 


w 7 0.l±i;S 


[-0.3,0.5] 


3.6 


4.6 O.Otg;? 



TABLE I: The signal yields and errors obtained from the sig- 
nal extraction fit, the ranges of observed biases from BB back- 
grounds, selection efficiencies (e), and the inferred branching 
fractions (B) for B° — > p°7, B + — * p + 7, and B° — > tu-y in the 
on-resonance data sample. The "Upper Lim." is a 90% C.L. 
limit. The efficiencies include the partial branching fractions 
for the p/uj decays considered. 



For the signal extraction, we perform an unbinned ex- 
tended maximum likelihood fit to the selected events. 
For B — > p7, the fit uses tties, AE* , and m„, whereas 
for B° — > wy, only m^s an d AE* are used. The mea- 
sured variables are largely uncorrelated, even after the 
p* w (or pl+^-^a) cut, allowing the probability density 
function (PDF) to be constructed as a product of in- 
dependent distributions for each variable. Since the BB 
backgrounds have PDFs that largely resemble continuum 
but are much smaller, the signal extraction uses only a 
continuum component to describe the background. Bi- 
ases due to BB backgrounds are considered below. The 
signal ibes and AE* distributions are described by the 
Crystal Ball shape with the exception of the ibes 
distribution for B° — > p°7, where the Gaussian distribu- 
tion is used. The relativistic Breit-Wigner lineshape is 
used for the signal distribution. The signal PDF 
parameters are obtained from MC simulation. The back- 
ground mES and AE* distributions are described by the 
ARGUS threshold function |l5j and a second-order poly- 
nomial, respectively. The background m v „ function is a 
sum of a Breit-Wigner component and a combinatorial 
component described by a first order polynomial. The 
background PDF parameters are determined in the fit, 
with the exception of the m^^ resonant fraction, which 
is fixed to the value measured in off-resonance data. 

The AE* vs. mss distributions of the selected B — > pry 
and B° — > uj^ candidates are shown in Fig. and the fit- 
ted signal yields are shown in Table [I] No significant sig- 
nal is seen in any mode. The quality of the fit is checked 
by comparing the overall likelihood of the fit with values 
obtained from an ensemble of parameterized MC simula- 
tions and found to be within the range expected. 

We consider three sources of systematic uncertainty in 
this analysis: the modeling of BB backgrounds, the sig- 
nal reconstruction efficiency, and the fixed parameters of 
the PDFs used in the fit. The first of these is "additive" 
in that it could result in background adding to the fit- 
ted signal yields. The last two are "multiplicative" in 
that they affect the way a given signal is interpreted as 
a branching fraction. 



The effect that BB backgrounds have on the fitted sig- 
nal yields is studied in parameterized MC simulations in 
which the BB background shape in the m ES -AE* plane 
is modeled with both one- and two-dimensional distribu- 
tions. Also, the rates of the dominant background modes 
are varied within wide ranges. For b — > sj (including 
B — > K*j), the normalization is varied between zero and 
twice the nominal value to conservatively account for un- 
certainties in kaon misidentification. For B + — » p + 7r° 
decays the branching fraction is varied between zero and 
twice the expected rate of 2 x 1CU 5 ^(|. Much lower 
branching fractions are expected for B° — > p°7r° and 
B° — » lutt [lif. so these cause negligible backgrounds. 
The small biases shown in Table [I] confirm that the BB 
PDFs are similar to those of continuum background. 

All signal-efficiency systematic uncertainties, except 
those related to the neural network and the uj mass, which 
are described above, are estimated in Ref. 0- The largest 
uncertainties, which arise from the neural net efficiencies, 
are 5%, 5%, and 10% for B° -> p°j,B+ -► p+ 7 , and 
B° — » 077 respectively. The ir° efficiency also contributes 
a 5% uncertainty to B + — > p + j and B° —> uj. 

The fixed parameters of the signal PDFs are stud- 
ied in fits to data for the topologically and kinemati- 
cally similar, but much more common, B — > K*j de- 
cays: B° -> *T*°7, K*° -> K+7T- for B° p° 7 
and B + -> K* + j, K*+ -> K + tt° for B+ -> p+ 7 and 
B° — » ujj. In these fits, the signal PDF parameters are 
allowed to float. The signal event yields are compared to 
those expected from the branching fractions measured in 
Ref. || and found to agree. 

The statistical uncertainties of the PDF parameters, 
one of which is the resonant fraction in the background 
TOtttt distribution, are used as ranges within which we 
vary the parameters of the B — > (p/uj)j fits. The result- 
ing variations in the fitted signal yield, which amount to 
5% for B° -> p°7 and B° -> w 7 and 10% for B+ -> p+7, 
are taken as systematic uncertainties. The total multi- 
plicative systematic error, including the signal efficiency 
uncertainty, is 8% for B° — ► p°j and 13% for B + — > p + 7 
and B° — > ujj. 

We assume S(T(45) -> B°B°) = B{T(AS) -» B+B~) = 
0.5. In calculating upper limits, we correct for bias from 
BB backgrounds by subtracting the smallest observed 
bias, which is found to be negative for all three modes, 
from the signal yield. We include the effects of the multi- 
plicative systematic uncertainties b y u sing an extension 
|17j of the method described in Ref. |lj| , wherein the sys- 
tematic and statistical errors are convolved. The result- 
ing 90% C.L. upper limits for the branching fractions are 
B(B° -> p° 7 ) < 1.2 x IO" 6 , B(B+ -> p+7) < 2.1 x 10~ 6 , 
and B(B° -> wy) < 1.0 x 10~ 6 . Although no signifi- 
cant signals are seen, Table U shows the measured B for 
each mode. For this calculation, we subtract a bias cor- 
responding to the center of the allowed range, treat the 
half-width of the range as the systematic error, and add 
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systematic and statistical errors in quadrature. 

We also calculate a combined limit for the generic pro- 
cess B — > p7 by assuming T(B — ► pj) = T(B + — > 
p + 7) = 2 x r(B° — » p°7) and using the lifetime ratio 
t b +/t b o = 1.083 ±0.017 9]. The resulting 90% C.L. up- 
per limit is B(B — » frj) < 1.9 x 10~ 6 . Using the measured 
value of B(B — > if*7) [8j, this corresponds to a limit of 
B(B -► pj)/B(B -> if* 7 ) < 0.047. 

This limit may be used to constrain the ratio of CKM 
elements |Vf<j/^t«| by means of the equation 0: 



B(£ -> if * 7 ) 



u, 



1/, 




C 2 [l + Ai?], 



where £ describes the flavor-SU(3) breaking between p 
and if*, and AR accounts for annihilation diagrams. Ai? 
is different for p° and p + , but we do not take this into 
account here. Both £ and Ai? must be taken from theory 
and there are several different 0, ^| values published. 
As an example, we choose the values £ = 0.76 ±0.10 and 
Ai? = 0.0 ± 0.2. We adjust both parameters down by 
one a and find the limit \V td /V ts \ < 0.34 at 90% C.L. 

In conclusion, we have found no evidence for the ex- 
clusive b —> dry transitions B —> pj and B — > oj 7 in 
84.4 ± 0.9 million BB decays studied with the BABAR 
detector. The 90% C.L. upper limits on the branching 
fractions are significantly lower than previous values and 
start to restrict the range indicated by SM predictions 
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